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^ ! Abstract 

O ' Stars and black holes are sources of gravitational radiation in many 



phases of their life, and the signals they emit exhibit features that are 
characteristic of the generating process. Emitted since the beginning 
of star formation, these signals also contribute to create a stochastic 
background of gravitational waves. We shall show how the spectral 
£NJ . properties of this background can be estimated in terms of the energy 

spectrum of each single event and of the star formation rate history, 
' which is now deducible from astronomical observations. We shall fur- 

O ■ ther discuss the process of scattering of masses by stars and black 

holes, showing that, unlike black holes, stars emit signals that carry a 
' clear signature of the nature of the source. 

as : 

o 

1 Introduction 

The extraordinary experimental effort done in recent years to design and 
| construct sophisticated and sensitive gravitational antennas, allows to hope 

^ ■ that gravitational waves will soon be detected. When this will happen, will 

we be able to identify the nature of the emitting source, or to decide whether 
a detected stochastic background was emitted during the early stages of the 
life of our Universe, or later, by a cosmological population of astrophysical 
sources? In order to answer these questions, waveforms and energy spectra of 
the signals emitted in different processess have to be computed and analysed 
to extract the specific information they provide on the generating process. 
In this paper we shall discuss some astrophysical phenomena in which this 
characterization is possible. 

Many phases of the life of a star are accompanied by the emission of grav- 
itational waves. For instance, a compact, rotating star emits radiation if it 
has a time varying quadrupole moment due to a triaxial shape or to a preces- 
sion of its angular velocity around the simmetry axes. When stellar objects 
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interact, the wave emission can be associated to the orbital motion, or to the 
excitation of their proper modes of vibration. Further on, intense bursts of 
radiation can be emitted in catastrophic events such as the gravitational col- 
lapse or the coalescence of binary systems. All these processes occurred since 
the beginning of star formation, and therefore they also contribute to create 
a stochastic background of gravitational radiation. To evaluate the spectral 
properties of the contributions of different sources we need a model of the 
energy spectrum emitted in each single event, as well as the rate of events of 
the selected type occurred in the recent or far past, i.e. the star formation 
rate history. As we shall later see, this important piece of information can 
be deduced from recent astronomical observations. 

Among the processes that are associated to the emission of gravitational 
waves, we shall first consider the final stages of the life of a sufficiently massive 
star which collapses to form a black hole. In this case numerical simulations 
show that, unless the collapse is dominated by very strong pressure gradients 
or by very slow bounces, the newborn black hole wildly oscillates and radi- 
ates its residual mechanical energy in gravitational waves, at well defined 
frequencies and damping times that are associated to its proper modes of 
vibration: the quasi-normal modes. We will then show how the energy spec- 
trum emitted in each event can be convoluted with the observation-based 
star formation rate history to compute the spectral energy-density of the 
resulting background, and we will show that it maintains memory of the 
quasi-normal mode signature of each collapse. This procedure is applicable 
to any cosmological population of astrophysical sources for which a model of 
the emitted energy spectrum is available. 

Another interesting astrophysical process that is associated to a well char- 
acterized signal is the scattering of a mass by the gravitational field of a mas- 
sive body, either a compact star or a black hole. We shall describe how the 
energy spectra and waveforms can be computed in the framework of a per- 
turbative approach, and show that the gravitational signals carry a specific 
information on the nature of the emitting source. 

2 The gravitational collapse to a black hole 

In 1939, Oppenheimer & Snyder QXJ studied the collapse of a spherical cloud 
of a pressureless gas of particles (dust), and showed that once the collaps- 
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ing cloud reaches the horizon, it cuts the communications with the external 
world and forms a black hole. In that case, due to the assumed spherical 
symmetry, no gravitational radiation emerged. However, gravitational waves 
are emitted if the collapsing configuration deviates from sphericity, and this 
is the basic idea of a perturbative approach introduced by C.T. Cunningham, 
R.H. Price and V. Moncrief 0, [|] and subsequently applied by E. Seidel and 
T. Moore |5], || (for a recent review see [0]). The perturbative approach con- 
sists in the following. Due to a sudden change in some equilibrium variable, 
a sperically symmetric compact star with an assigned equation of state col- 
lapses. Typically, the pressure is changed so that the star is no longer in 
equilibrium. The collapsing configuration is obtained by numerically inte- 
grating the Einstein equations coupled to the equations of hydrodynamics, 
and this exact solution - 'exact' in the sense that it is a solution of the fully 
non-linear equations - is used as a background about which linear perturba- 
tions are considered. The perturbed equations, derived in the framework of a 
gauge invariant formalism |2j, |J , are then numerically integrated to compute 
the energy spectrum and the waveforms. 

Alternatively, the gravitational collapse can be studied by directly inte- 
grating the fully non-linear equations of gravity, and at present there is only 
one fully relativistic numerical simulation of collapse to a black hole. In 1985 
R. Stark and T. Piran |J numerically integrated the equations describing the 
collapse of a rigidly rotating, axisymmetric polytropic star. Although this 
simulation is based on some simplifying assumptions, it can be considered as 
a good model to extract information on the collapse of a massive star to a 
black hole; indeed, it has be shown that the spectrum of the emitted grav- 
itational energy exhibits some distinctive features that are, to some extent, 
independent of the initial conditions and of the equation of state of the col- 
lapsing star, and that depend only on the black hole mass and on its angular 
momentum. 

The dynamics of the collapse proceeds as follows. If the total angular 
momentum a = is greater than a critical value a cr u the rotational 

c 

energy dominates, the star bounces and no collapse occurs. Conversely, if 
a < a cr u a black hole forms. In figure la we plot the waveform, i.e. the 
/i + -component of the metric tensor of the gravitational wave emerging at 
infinity, evaluated in the TT-gauge. The h x -component is much smaller 
and does not significantly contribute to the emitted radiation. The function 
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Figure 1: The h + component of the gravitational radiation (a) and the 
average energy flux per unit frequency (b ) emitted during the axisymmetric 
collapse of a rotating, polytropic star to a black hole, are plotted for assigned 
value of the angular momentum j^j. These data refer to a collapse ignited 
by a pressure reduction to a fraction f p = 0.01 of the equilibrium central 
pressure. 
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f{v) plotted in figure lb is an average energy flux computed as follows 

The amplitude of f(y) scales as the fourth power on the angular parameter 
a, and it is proportional to the square of the mass of the collapsing core; the 
efficiency of the process is AE GW / M core c 2 < 7 x 10~ 4 . 

Figure lb shows that if the angular momentum increases more energy 
is emitted in the low frequency region. This is due to the fact that for 
high values of a the star becomes flattened into the equatorial plane, and 
bounces vertically before collapsing. However, the more distinctive feature 
of these energy spectra is the peak which occurs at a frequency that depends 
on the angular parameter; for instance v max — 8-7 for a = 0.79, and 
v max ~ 9.0 kHz for a = 0.94. Moreover, the ringing tails exhibited by 
the waveforms plotted in figure la are damped sinusoids at these frequencies. 
In order to understand what is the physical phenomenon which associated 
to the emission of these peaks, we need to define the quasi-normal modes of 
black holes, that emerge from the theory of black hole perturbations 0. 



2.1 The perturbations of a rotating black hole 



It is known that the radial evolution of the perturbations of a rotating 
black hole is described by a single master equation for a radial wavefunc- 
tion Ri m (r,u): 



A 



d 2 R 



Or 



lm + 2(s + l)(r - M)^ + V(r, u)R lm = 0, 



(2) 



where 

V(r,u) 



A' 



- Mu{r 2 - a 2 
A = r 2 - 2Mr + a 2 



(r 2 + a 2 ) 2 uj 2 - AaMrmuj + a 2 m 2 + 2is amir - M) (3) 
2iscor — a 2 oj 2 — Ai. 



+ 



Hm 



s is the spin- weight parameter, s = 0, ±1, ±2 respectively for scalar, elec- 
tromagnetic and gravitational perturbations, and the sign indicates whether 
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the ingoing or outgoing radiative part of the field are considered. Ai m is a 
separation constant. This wave equation was derived by S. Teukolsky in 1972 



[ PH| , |TT| ]. The quasi-normal modes of a black hole are defined to be complex 
frequency solutions of the wave equation |2|, that behave as a pure outgoing 
wave at infinity, where radiation has to emerge, and as a pure ingoing wave 
at the horizon, since nothing can escape from a black hole horizon. The 
eigenfrequency of a mode is complex because the emitted radiation damps 
the oscillations, and the real part is the pulsation frequency, whereas the 
imaginary part is the inverse of the damping time. Numerical integrations 
of the wave equations have shown that these frequencies are characteristic of 
many processes involving dynamical perturbations of black holes, and that 
an initial perturbation decays, during its very last stag superposition 
of these pure modes. In table 1 we tabulate the characteristic frequency 
of the lowest quasi-normal mode of a rotating black hole, for I = 2 and 



m = 0, computed by Leaver |12"||. The corresponding values in physical units 



can be evaluated by using the following formulae 

c(Mu ) 32.26, , 7TT 
27m ■ M Q n 
nM Q _ n ■ 0.4937 • KT 5 
{MuJi)c ~ (Mui) S ' 

where M Q = 1.48 ■ 10 5 cm and we have assumed that the black hole mass 
is M — nM Q . Thus we easily find that for I = 2 

a = 0.79 vq = 8.6 kHz, a = 0.94 u = 9.0 kHz. 

If we now go back to the energy spectra plotted figure lb, we see that they 
peak at a frequency which is very close to uq, i.e. to the frequency of the 
lowest m = quasi-normal mode of the rotating black hole which forms 
as a result of the collapse. It should be mentioned that the excitation of the 
quasi-normal modes also emerges in simulations of the collapse to a black hole 
done in the framework of perturbative approaches. Thus, the quasi-normal 
mode peak is a specific signature of the emitted energy spectrum, since it 
indicates that the newborn black hole oscillates in these modes, releasing its 
mechanical energy in the form of gravitational waves. 

We may ask to what extent is this signature a general feature of the 
signals emitted during the collapse of a sufficiently massive star to a black 



6 



Table 1: The characteristic frequency of the first quasi-normal mode of a 
rotating black hole is tabulated, for 1 = 1 and m = 0, for different values 
of the angular parameter a. 



a 


Mu + iMuJi 


0.0 


0.3737+i0.0890 


0.2 


0.3751+i0.0887 


0.4 


0.3797+i0.0878 


0.6 


0.3881+i0.0860 


0.8 


0.4019+i0.0822 


0.9 


0.4120+i0.0785 


0.98 


0.4223+i0.0735 


0.9998 


0.4251+i0.0718 



hole. A study done by E.Seidel using the perturbative approach fl3|l , has 
shown that if the collapse is dominated by extremely high pressures or by 
slow bounces, the infalling material can be considerably slowed down and the 
formed black hole may be left with very little residual energy to radiate away 
in gravitational waves. In this case the quasi-normal modes peak may not be 
as pronounced as in figure lb. However, these are quite extreme situations, 
and it is reasonable to presume that the gravitational collapse will leave 
behind a newborn black hole with some residual mechanical energy, which 
will be radiated in gravitational waves at the frequencies of the black hole 
quasi-normal modes. The efficiency of the gravitational emission evaluated 
by Stark and Piran has to be taken as an indication: it may be smaller, if 
strong pressures and slow bounces dominate the collapse, but it may also be 
higher, if the collapse is non-axysimmetric. 



3 From a single event to a stochastic back- 
ground 

The gravitational signals emitted by newly formed black holes throughout 
the Universe superimpose to form a stochastic background. In order to eval- 
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uate the spectral properties of this background and verify whether it keeps 
memory of the generating process, besides the energy spectrum emitted in 
each single event, we need to know the rate of events as a function of the 
cosmological redshift. This information can be deduced from the observative 
data collected by the Hubble Space Telscope, Keck and other large telescopes 
[|P|j-||15|, which, together with the completion of several large redshift surveys 
[|L6[|-||17| have enabled to derive coherent models for the star formation rate 
evolution up to redshifts of ~ 4 — 5. 

In figure || we plot a recently proposed fit of the data obtained at different 
redshifts, representing the star formation rate density p\{z), which is the 
mass of gas which forms stars per unit time and comoving volume element 
|]r5). The observative data-points use the rest-frame UV-optical lumi- 



nosity as an indicator of the star formation activity in distant galaxies, and 
it is known that even a relatively small amount of dust can absorbe the UV 
light and reradiate it in the far-IR. This causes an attenuation of the ob- 
served UV-luminosity (dust extinction) and leads to an underestimate of the 
real star formation activity. To account for this problem, the data shown 
in fig. |H have been corrected according to the Calzetti dust extinction law 
(see [0, \Wjj). The amount of dust correction to be applied at intermediate- 
to-high redshift 2-1 — 2, is still relatively uncertain. However, since the 
energy emitted by a single source decreases as the square of the inverse of 
the luminosity distance, the gravitational wave backgrounds generated by 
astrophysical sources are almost insensitive to the high redshift behaviour of 
the star formation rate. Using the function p\(z) plotted in fig. |3|, the rate 
of core-collapse supernovae can be computed as follows |HJ 



R ^ z ) = r dz 'b^l^ [ dM'^(M'). (5) 

v ' Jo 1 + z' dz Jam v ' v ' 



P*{z') dV 
1 + z' 

The the factor (1 + z)~ x takes into account the dilution due to cosmic ex- 
pansion and $(M) is the initial mass function chosen to be of Salpeter type, 
$(M) oc M~( 1+x *> with x = 1.7. 

The mass range, AM, depends on the nature of the considered source. Nu- 
merical studies on stellar evolution have shown that single stars with mass 
exceeding 8M evolve through all phases of nuclear burning, ending their 
life as core collapse supernovae (this class includes Type II and Types Ibc 
supernovae). While there seems to be a general agreement that progenitors 
with masses in the range 8M Q ~ M ~ 20M Q leave neutron star remnants, the 



8 



bo 
o 




Figure 2: The star formation rate density in units of M Q yr _1 Mpc -3 is plotted 
as a function of redshift for a cosmological background model with Qm = 1, 
£Ia = 0, Hq = 50kms _1 Mpc _1 and a Salpeter initial mass function. 
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value of the minimum progenitor mass which leads to a black hole is still 
uncertain, mainly because of the unknown amount of fall back of material 
during the supernova explosion [21], [22]. In the following we shall assume 
that the lower threshold for black hole formation is M~25 M Q . 

The rate of events ([5]) depends also on the cosmological background 
through the comoving volume element 

dV = 4vr (±) r 2 1 dz, (6) 

\H J . s J(l + z) 2 (l + tt M z)-z(2 + z)n A 

where r is the comoving distance 

c 



S(^K J\l + z') 2 {l + n M z')-z'{2 + z')n A \- 1 l 2 dz'), (7) 



and the function S is given by |2 



Vl m + Vl a >1 S(x)=am(x) Q K = 1 - Q AI + tt A (8) 
M + < 1 S(x) = sinh(a;) Q K = Q M + Q A — 1, 
VL M + Q A = 1 S(x) = x Q K = 1. 

If we consider the following cosmological backgrounds [] 

A) tt M = l tt A = h = 0.5 (9) 

B) tt M = 0.3 Q A = 0.7 h = 0.6 

C) tt M = 0.4 tt A = = 0.6 

the total rate of core-collapse supernovae leading to a black hole or to a 
neutron star ranges, respectively, between 

Rbh ~ 3.3 — 4.7 events/ s, Rns ~ 13.6 — 19.3 events/ s. (10) 

The main difference between the three cosmologies is introduced by the geo- 
metrical effect of the comoving volume, and it is significant at 2-1 — 2. This 



1 We shall assume 

h = — — km/s/Mpc, 
100 ' ' y ' 

where H is the Hubble constant. 
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implies that the gravitational backgrounds, which are mainly contributed by 
sources at z~l — 2, are almost insensitive to the cosmological parameters. 
It is now possible to compute the spectral energy density of the gravitational 
wave background produced by an assigned population of sources 



dE (y) 
dtdSdv 



oo JAM 



f{v)-dR{z,M) 



(11) 



where dR(z,M) is the differential rate, and f{y) is the average energy 
flux per unit frequency, chosen as a model for the class of sources under 
consideration, emitted by a source located at a luminosity distance di,{z) = 
(1 + z) r, 

m = nV« r ^ f ( d T^r) siMd - ( 12 ) 

Av:d L (z) 2 Jo Jo \dvdVt ) 

For instance, for core collapses to a black hole we shall use as a model of f(v) 
the energy spectrum shown in figure lb, rescaled with the luminosity 
distance and suitably redshifted. 

From eq. flll]) we can further derive the closure energy density of gravi- 
tational waves per logarithmic unit frequency, 



n 



dp, 



GW 



GW 



V 

C 3 p c 



dE, 



GW 



(13) 



where p c 



3HI 
SttG 



S h (v) 



1 

p c dlogv c 3 p c dtdSdv 
1.9 x 10 -29 h 2 g cm~ 3 , and the spectral strain amplitude 



2G 

7TC 3 ^ 2 



dE, 



GW 



dtdSdv 



1/2 



f 3 H 
M7T 2 Z/ 3 



■ n 



GW 



1/2 



(14) 



which is the quantity to be compared with the detectors sensitivity. The 
spectral properties of the background produced by a cosmological populations 
of core-collapse supernovae leaving behind a black hole, computed by this 
procedure, are shown in figure [I]. We plot the closure energy density and the 
strain amplitude for three selected values of the angular parameter of the 
formed black holes, assumed to rotate all at the same speed since we do not 
know the distribution of black holes angular momenta. 

The strain amplitude exhibits a sharp peak at a frequency which ranges 
within ~ [230 — 340] Hz, depending on the chosen value of the angular mo- 
mentum, and with a maximum amplitude y/Sh max ~ [1 — 8] x 10~ 28 Hz^ 1 ^ 2 . 



11 



500 1000 1500 2000 2500 3000 3500 4000 4500 5000 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 

v (Hz) v (Hz) 



Figure 3: The closure energy density Slew (left) an d the spectral strain 
amplitude S 1 ^ 2 (right) generated by a cosmological population of core col- 
lapse supernovae leaving behind a black hole, is plotted as a function of the 
observational frequency, for three values of the angular parameter a and for 
a flat cosmology with vanishing cosmological constant. 
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This peak is reminiscent of the excitation of the quasi-normal modes of the 
formed balck holes, showing that the spectral properties of this background 
keeps memory of the generating process. 

It should be mentioned that since according to eq. ( |10D the event rate is 
of the order of 3 — 5 events/ s, and since a signal emitted in a collapse to 
a black hole has a very short durations (cfr. eq. typically of the order of 
a few milliseconds, the background generated by newborn black holes has a 
shot noise character. 

The procedure to determine the characteristics of a stochastic background 
described in this section can be applied to any population of astrophysical 
sources for which a model for the energy spectrum is available. For example, 
in ref. we have applied this method to study the background produced 
by a cosmological population of young, rapidly rotating neutron stars that, 
due to the r-modes instability, are expected to radiate a large fraction of 
their rotational energy in gravitational waves. A gravitational background 
can be detected by cross-correlating the output of two gravitational antennas 
for a sufficiently long time interval (typically one year), and in refs. |2(| and 
24|] we have computed how interferometric and resonant detectors operating 



in coincidence might respond to our predicted astrophysical backgrounds. 
We find that whereas the sensitivity of the first generation of interferometric 
antennas like VIRGO and LIGO will be too low to detect these signals, the 
planned sensitivity of the advanced version of these experiments would allow 
the detection, provided two interferometers with similar characteristics are 
located nearby. 



4 Gravitational waves from stars: a pertur- 
bative approach 

A number of numerical simulations of the gravitational collapse to a neutron 
star have shown that, unlike the case of collapses to a black hole, the gravi- 
tational signals computed for these processes strongly depend on the initial 
conditions, on the equation of state of the collapsing star, and on the de- 
tails of the collapse P5[]-P7|; for this reason, a model of the energy spectrum 
emitted in a collapse to a neutron star is still not available. The situation 
changes if we consider processes occurring after a neutron star has formed, 
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and we shall now show how some interesting information can be derived by 
using a perturbative approach. 

The theory of stellar perturbations was formulated in the framework of 



General Relativity by Thorne and his collaborators since 1967 |[28||-||33||, and 
it was successfully applied to determine the frequencies of the quasi-normal 
modes of oscillation for a wide range of stellar models [|34], ^] . Much recently, 
the theory has been reformulated in close analogy to the theory of black hole 
perturbations 36||-[43[|, and new phenomena have emerged, which do not 



have a newtonian counterpart. A detailed discussion of the many interesting 
aspects of this theory is beyond the scope of this paper, where we want to 
focus essentially on the characteristics of the gravitational signals emitted 
in astrophysical processes. Thus, we shall rather consider an application of 
the theory to a specific process, and show that the emitted signals exhibit a 
clear signature of the nature of the source. In particular, we shall consider 
a mass m which, interacting with the gravitational field of a star of 
mass M deviates from its original trajectory moving toward the star, 
reaches a periastron (the turning point) and then moves away in an open 
orbit. Under the assumption that itlq << M, the problem can be solved 
by using a perturbative approach. The matching of the solutions of the 
perturbed equations inside and outside the star constitutes the most delicate 
technical point MM . Indeed, while the interior solution is found by integrating 



the equations for the perturbed metric tensor coupled with the equations of 
hydrodynamics, the exterior solution cannot be found by using the same 
tensorial approach because the source term of the equations for the scattered 
mass diverges at the turning point. However this difficulty can be overcome 
by introducing a suitable wavefunction, related to the Weyl scalar ^4, 



which carries information on the radiative part of the field. In ref. |45[ the 
equations describing the perturbations of a star induced by a scattered mass, 
have been integrated for a star with a polytropic equation of state p = Kp n , 
with K = 100 km, n = 2, and central density p c = 3 ■ 10 15 g/cm 3 . 
The radius and mass of the star are, respectively, R = 8.86 km, and 
M = 1.266 M©, with a ratio R/M = 4.7 The energy spectrum emitted 
in gravitational waves has been computed for the following sets of orbital 
paramenters 

• a) L = 5 and E = 1.007, which correspond to a turning point 
located at r t = 9.2 M 
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Figure 4: The £ = m = 2 - component of the energy spectrum of the 
gravitational radiation emitted when a mass m is scattered by a polytropic 
star of mass M, is plotted versus the normalized frequency cuM . The 
orbital parameters L and 2£ are chosen in such a way that the turning 
point is at r t = 9.2 M (a) and r t = 5 M (b) 
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• b) L = 5 and E = 1.097, so that the mass can get closer to the 
star and r t = 5.0 M, 

and it is plotted in figure f| versus the normalized frequency u M. Q L = 
L z /M is the angular momentum of the scattered mass m , normalized to 
M, and E is its energy per unit mass. 

Figure [| shows that the energy spectra exhibit well defined peaks located 
at some particular frequencies, and in order to understand the physical pro- 
cesses that underlie this structure we need to mention what are the character- 
istic frequencies at which a star can oscillate and possibly emit gravitational 
waves. 

In newtonian theory, the classification of the modes of oscillations of a 
star is based on the behaviour of the perturbed fluid. The hydrodynamical 
equations show that when a star is perturbed each element of fluid moves 
under the competing action of two restoring forces, one due to the eulerian 
change in the density 5p, the other due to a change in pressure 5p. The 
modes are classified accordingly: g modes, if the prevailing driving force is 
due to Sp, p modes, if it is due to 5p. The two classes of modes occupy well 
defined regions of the spectrum, and they are separated by the fundamental 
mode, the f mode, characterized by having an eigenfunction that has no nodes 
inside the star. In General Relativity stellar oscillations are damped by the 
emission of gravitational waves, and in addition to the fluid modes there exist 
modes of the radiative field. Indeed, it can be shown that fluid motion is 
either negligible or totally absent at the corresponding frequencies. They are 
named w modes, characterized by high frequencies and short damping times 
HH , and s modes, that are slowly damped, do not excite any motion in the 
fluid and are characteristics of ultracompact stars (R/M-2.Q) p9[ . 



Thus, stars are characterized by a very rich set of possible modes of 
vibration, and it is interesting to ask whether these modes can be excited in 
real astrophysical processes, and how much energy they carry. For example, 
for the model of star considered in figure || the frequency of the fundamental 
mode is to M = 0.120, the first p -modes is u M = 0.266 and the first 
w-mode is u M = 0.53. The energy spectra of figure |] exhibit a pronounced 
peak at the frequency of the f mode, and more peaks appear, corresponding 



2 In figure || we plot only the £ — m — 2 component of the energy spectrum, because 
it provides the dominant part of the radiated energy. 
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(15) 



to the excitation of the first p -modes, if the scattered mass is allowed to get 
closer. 

It is interesting to compare this behaviour with that of a black hole per- 
turbed by a small mass in a similar scattering process (see j|7| for an ex- 
tensive review). It turns out that the black hole is rather insensitive to 
these processes: the energy is emitted essentially by the scattered mass as 
a synchrotron radiation, and most of it is radiated when the mass transits 
through the turning point. Indeed, the energy spectrum is peaked at a fre- 
quency which is related to the angular velocity of the mass at the turning 
point as follows 

K at ' r=r t 

In our case a) (r t = 9.2M), the frequency corresponding to the angular 
velocity of the mass at the periastron for 1 = 1 is u rt M = 0.092. The 
spectrum shown in figure |]a has a peak at that frequency, showing that part 
of the energy is still emitted by the mass as a synchrotron radiation, but 
the the very sharp peak which occurs at the frequency of the fundamental 
mode is dominant. A similar situation arises when the mass gets closer 
to the star (case b, r t = 5M, u> n M = 0.221), but again the emission is 
strongly dominated by the excitation of the modes of the star. Thus, we 
can conclude that, unlike black holes, whose quasi-normal modes are not 
excited in scattering processes, a mass moving around a star in an open 
orbit can excite the fluid modes of the star, to an extent that depends on 
how close is the encounter. The w-modes do not appear to be significantly 
excited in these processes, and since the contribution of the fundamental 
mode appears to be the dominant one, the emitted gravitational signal is a 
pure note corresponding to that frequency. 



5 Concluding Remarks 

The study of the gravitational radiation emitted by astrophysical sources 
is an open field to explore, and it requires the cooperative effort of scien- 
tists working in different fields. Complex phenomena like the gravitational 
collapse of a massive body or the coalescence of binary systems are paradig- 
matic: powerful computers and specific expertise in numerical techniques to 
deal with the strong field regimes typical of these processes are needed, as 
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well as experience in the physics of phenomena occurring at supernuclear 
densities and a deep mathematical knowledge of the equations of gravity. In 
addition, since to detect a signal it is important to know not only the wave- 
form and the energy it carries, but also the rate of occurrence, the knowledge 
of the sources distribution and evolution throughout the Universe, either the- 
oretical or based on observations, is also crucial. 

Although most of the interesting phenomena that are associated to the 
emission of gravitational radiation are highly non linear, much can be learnt 
by the use of approximation schemes. For instance the perturbative approach 
proves extremely useful in providing a physical understanding of many pro- 
cesses; it helps to interpret and confirm the results of fully non linear simu- 
lations, and to identify physical mechanisms operating in different regimes. 
An example of this synergy of the approaches is given in section 2, where the 
peaks of the energy spectrum emitted in a collapse to a black hole, computed 
by a fully relativistic numerical simulation, have been interpreted in terms 
of the excitation of the quasi normal modes predicted by the theory of black 
hole perturbations. 

Furthermore, we have shown how theory and observations can be matched 
together to evaluate the spectral properties of the background of gravitational 
waves produced by these sources: the energy spectrum of each single event 
derived by numerical simulations has been convoluted with the star formation 
rate history deduced from astronomical observations. 

The work presented in this paper can be extended in several directions. 
We plan to apply the procedure developed to evaluate the background of 
gravitational waves to other cosmological populations of astrophysical sources, 
like binary systems formed by neutron stars, black holes and white dwarfs 
48fl . Furthermore, we shall compute energy and waveforms of the signals 



emitted by stars perturbed by masses orbiting in closed orbits, and investi- 
gate the relation between the amplitude of the emitted wave and the equation 
of state prevailing in the star interior, for different star models ranging from 
sun-like stars to white dwarfs and neutron stars. Finally, we plan to investi- 
gate the role that the excitation of the quasi-normal modes can play in the 
coalescence of binary systems, and extend all these calculations to rotating 
stars. 
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